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Abstract: We revisit the strong coupling limit of the cusp anomalous dimension in 
planar J\f = 4 super Yang-Mills theory. It is known that the strong coupling expansion is 
asymptotic and non-Borel summable. As a consequence, the cusp anomalous dimension 
receives non-perturbative corrections, and the complete strong coupling expansion should 
be a resurgent transseries. We reveal that the perturbative and non-perturbative parts 
in the transseries are closely interrelated. Solving the Beisert-Eden-Staudacher equation 
systematically, we analyze in detail the large order behavior in the strong coupling pertur¬ 
bative expansion and show that the non-perturbative information is indeed encoded there. 
An ambiguity of (lateral) Borel resummations of the perturbative expansion is precisely 
canceled by the contributions from the non-perturbative sectors, and the final result is real 
and unambiguous. 
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1 Introduction 

In recent years, there is remarkable progress in the understanding of the AdS/CFT duality 

[1-3]. In planar AdS/CFT between AT = 4 super Yang-Mills (SYM) theory and string 
theory on AdS§ x S 1 * * * 5 , integrability is a key concept in the investigation. Here we refer 
to [4] for a comprehensive review on the AdS/CFT integrability, but of course there are 
many important results after that review. In particular, recently a new formulation, called 

Quantum Spectral Curve, to solve the spectral problem in planar AT = 4 SYM was proposed 
in [5, 6], 

In this work, we focus on a well-studied quantity called the cusp anomalous dimension 
T cusp (g) (or sometimes called the scaling function). Throughout this paper, we use the 
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conventional notation of the coupling constant g , which is related to the ’t Hooft coupling 

A = gy M N b y 

VA 


g = 


47T 


( 1 . 1 ) 


The cusp anomalous dimension appears in several contexts. It appears as a UV divergence 
of a (light-like) Wilson operator with a cusp or as an IR divergence of a gluon scattering 
amplitude. These divergences are closely related by the so-called Wilson loop/amplitude 
duality [7-10]. It is well-known that the cusp anomalous dimension also appears as a 
logarithmic divergence of conformal dimension A for a twist-two operator with large spin 
S: 

A - S = 2T cusp (g)logS-\ -, S -A oo. (1.2) 


Note that this logarithmic behavior is universal for any g, thus the cusp anomalous dimen¬ 
sion is a good interpolating function from weak to strong coupling. From these examples, 
it is obvious that to understand the cusp anomalous dimension is an important task. 

Very surprisingly, the cusp anomalous dimension in planar Af = 4 SYM can be com¬ 
puted at any coupling by solving the so-called Beisert-Eden-Staudacher (BES) equation 
[11] (see also [12])! This is one of the greatest achievements in the AdS/CFT integrability. 
We can learn many things through this equation. At weak coupling, the cusp anomalous 
dimension admits the standard perturbative expansion in g 2 : 


r c us P (sO — 


ttV 11t rV / 73-7T 6 

3 45 V 630 


4C(3) 2 ) g 6 + 0(g s ) 


g o. (1.3) 


As observed in [11], this weak coupling expansion is a convergent series with finite radius 
g = 1/4. Therefore, it is expected that there are no non-perturbative corrections 1 of 
the form e~ A ^ 9 , and the perturbative expansion (1.3) is sufficient to reconstruct the full 
function r cusp (g'). The weak coupling result was confirmed up to four loops [16, 17]. 

At strong coupling, the problem is much more involved. The BES equation predicts 
the strong coupling expansion 

/ \ ( 3 log 2 K \ 

r„p(9) = 2» (l - — -^+ -). I 1 ' 4 ) 

where K = 0.915965594... is Catalan’s constant. These coefficients were first predicted in 
[18] by the numerical analysis. In [19-21], the leading coefficient was computed analytically, 
and then, in [22], the expansion was systematically computed up to 1/g 40 . Of course, the 
strong coupling prediction (1.4) should be compared with the direct string worldsheet 
computation. In fact, the worldsheet computation up to two loops shows that the strong 
coupling prediction (1.4) is perfectly reproduced [23-26]. This agreement is obviously a 
strong evidence of the planar AdS/CFT duality. 

lr To be precise, a convergent series is not a sufficient condition for non-existence of non-perturbative 
corrections. A counter-example is the exact planar free energy in ABJM theory [13]. As shown in [14], it 
receives a non-perturbative correction of the form e ~ 27r '/ 2 >' a t strong coupling even though its perturbative 
1/VA expansion is convergent. The same also happens for the so-called interpolation function [15] (see 
(6.3)). 
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However, this is not the end of the story at strong coupling. A crucial observation in [22] 
is that the strong coupling perturbative expansion (1.4) is very likely asymptotic and non- 
Borel summable due to singularities on the positive real axis in Borel plane. At first glance, 
the non-Borel summability causes a serious problem in the perturbative resummation, and 
it seems that the Borel resummation procedure does not work any more. However, there is a 
beautiful resolution for this resummation problem. As well-known, the Borel singularities 
on the positive real axis cause an ambiguity of Borel resummations because one has to 
avoid the singularities when performing the inverse Borel transform (see figure 1). There 
are several choices for the integration contours to avoid the singularities. The important 
point is that the ambiguity has non-perturbative order e~ Aa , where A is related to the 
closest Borel singularity to the origin. Therefore a natural expectation is that there is a 
non-perturbative correction of the same order, and that its Borel resummation precisely 
cancels the ambiguity in the perturbative resummation. As a result, the total sum should 
give the same answer for any choice of the contours in the Borel resummations. 

This beautiful structure of cancellation of ambiguities goes under the name of resur¬ 
gence and it was developed by Ecalle [27]. Since then it has been applied at first in quantum 
mechanical systems [28, 29] and only very recently it has been applied to quantum field 
theory [30-32] to obtain a weak coupling interpretation of the IR renormalons . 

Similarly, in the cusp anomalous dimension, the Borel resummation of the strong cou¬ 
pling perturbative expansion alone is insufficient to reconstruct the full function T cusp ((7). 
To resolve the Borel resummation problem, we need non-perturbative corrections. In a se¬ 
ries of works [22, 33-35], it turned out that the cusp anomalous dimension should receives 
the following non-perturbative corrections 

1 =^ = r£ p ( 9 ) - JA rw p ( 9 ) + I A ? if p ) p ( 9 ) + o(A«), g -»(1.5) 

where Tcusp(fi') is the perturbative contribution above and A 2 is a non-perturbative scale, 
related to the’t Hooft coupling as follows 

A 2 = u^(27t (1.6) 

Note that the (complex) parameter a depends on a choice of Borel resummations, as 
will be seen later. The non-perturbative scale is closely related to the mass gap of the 
0(6) sigma model 2 , which describes the effective string worldsheet theory in the high spin 
limit, as explained in [36] (see also [37]). In each non-perturbative sector, Tcusp(g) has the 
asymptotic (non-Borel summable) 1/g expansion 3 : 

00 -p(n) 

■^) = E^c < L7 > 

2 The 0(6) sigma model, being an asymptotically free two-dimensional field theory, is affected by the so 
called IR renormalons. See [38] for a recent discussion of the resurgence properties of the O(N) models and 
the connections with the IR renormalons problem. 

3 One can always set rj)” 3 = 1 by factoring out an appropriate factor. 
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(n) 

In principle, one can compute the coefficients T) by solving the BES equation at strong 
coupling, but this is a highly non-trivial task. A very first few coefficients in the leading 
non-perturbative sector TcuspCg) were computed in [35]. 

One of our goals in this paper is to show that there are non-trivial relations among the 
coefficients rj"\ In particular, the perturbative coefficients rj 0 ^ already know information 
on the non-perturbative sector! See (5.2) for a concrete result. This is clearly a conse¬ 
quence of the resurgent analysis [39], recently developed in many contexts in theoretical 
physics. To confirm this fact, we develop the method in [35], and numerically compute the 
coefficients rj/^ (n = 0,1,2) up to £ = 180 with sufficiently high precision. The fact that 
the perturbative and non-perturbative parts are interrelated to each other provides us a 
strong consistency test of the strong coupling solution. If an obtained solution is wrong, 
one would encounter a discrepancy with the resurgent expectation. Moreover, one can pre¬ 
dict the non-perturbative correction from the perturbative result without solving the BES 
equation. We also demonstrate that the ambiguity of the Borel resummations in the per¬ 
turbative sector is precisely canceled by the same ambiguity in the non-perturbative sector. 
As a consequence, the final answer is always real-valued and unambiguous. Our results 
show that the cusp anomalous dimension at strong coupling is a resurgent transseries. 

The organization of this paper is as follows. In Section 2, we briefly review some 
basic concepts in Borel resummation and resurgence theory, we refer to [39] for a longer 
discussion. We discuss then, in Section 3, how to obtain the strong coupling expansion of 
the cusp anomalous dimension from the solution to the BES equation. To better understand 
the resurgent properties of the cusp anomaly in AT = 4, in Section 4 we analyze a toy model 
solution to the BES equation previously presented in [35]. In Section 5, we finally study the 
cusp anomaly in A/" = 4 and show how to reconstruct all the non-perturbative corrections, 
i.e. the full transseries solution for r cusp (( 7 ), simply from the perturbative strong coupling 
asymptotic expansion. We show how the transseries representation for T cusp (( 7 ) is free from 
ambiguities and we comment on the relation between the cusp anomaly and the mass gap 
of the 0(6) sigma model. We draw our conclusions in 6. 

Note added: While this work was being completed, we became aware of related work 
[40], which has some overlap with this paper. 

2 Borel resummations and resurgent transseries 

As we will review in details later on, in the strong coupling regime the perturbative con¬ 
tribution to the cusp anomaly, Tcu S p(g), and the perturbative corrections in each non- 
perturbative sector, Tcusp (g), all take the form of asymptotic power series in 1 /z (z = 2 ng) 

OO n 

/(*) = C + E]£l, z ^°°’ (2- 1 ) 

n=0 

where the coefficients f n diverge factorially like T(n + a), for some constant a, denoting a 
Gevrey-1 type series. For the terminology in the Borel analysis, see [39], for instance. 
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For this kind of asymptotic power series, the standard Borel transform 1 proves to be 
an extremely useful object, which is defined as 


= £ 


r(rc+ 1) 

n —0 


( 2 . 2 ) 


Note that this series converges at least in a neighborhood of the origin t = 0. We can thus 
obtain an analytic continuation of the original formal power series by a Laplace transform 

of B[f] 

/*e 10 oo 

Sof(z) = c+ d te~ tz , (2.3) 

Jo 

usually called Borel resummation of f(z ) in the direction 9. 

When B[f](t) contains singularities along a particular direction 9 , in the complex t- 
plane, also called Borel plane, we will say that 9 is a Stokes line for / (or equivalently 
B[f}). Along a Stokes line we cannot directly use the resummation Sq, but we can easily 
dodge the singularities by defining the two lateral resummations 


Sg+f{z) = c + 


Sg-f(z) =c + 



d te- tz B[f](t), 
dte- tz B[f}(t). 


(2.4) 


Note that if the direction 9 does not contain any singularity, the two lateral summations 
coincide with each others and with the standard Borel resummation (2.3): Sg = Sq+ = Sq- . 

When 9 is a Stokes line, we have Sq+ ^ Sq-, but it is still possible to relate the 
two different analytic continuations of the asymptotic series f(z ) via the so-called Stokes 
automorphism, 6q, in the direction 9 


S g + = S e - o6 e = S g - o (Id - Discs) , 

Sq+ — Sq- = —Sq- O Discs . 

Where Disc# encodes the full discontinuity across 9. 

By a simple contour deformation, we can rewrite the difference between the two re¬ 
summations along 9 + and 9~ as a sum over Hankel’s contours, and the discontinuity of S 
across 9 is given as an infinite sum of contribution coming from each one of the singular 
points, see Figure 1. 

From simple residue theory, it is easy to see that the difference between the two re- 
summations generates a new type of non-analytic term, beyond the realm of formal power 
series, schematically of the form 


(2.5) 


(Sq+ — Sq- ) ~ 2-7ri e ** * (a + 0(1/z)) , 


( 2 . 6 ) 


where the point t* is a singularity of B[f](t) in the direction 9, i.e., arg(i*) = 9. These 
exponentially suppressed terms are nothing new, in semi-classical calculations we know 

4 If the coefficients /„ grows asymptotically as T(n + a), it is sometimes more useful to use a slightly 
different Borel transform where the coefficient f n gets divided precisely by F(n + a). 
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Figure 1. The difference between left and right resummation along the singular direction 6 as a 
sum over Hankel contours. 


that generically the perturbative expansion receives non-perturbative corrections usually 
referred to as instantons corrections. So, whenever we have to deal with an asymptotic 
perturbative expansion, we expect that the most general solutions to our problem, in the 
case at hand the cusp anomalous dimension at strong coupling, will take the transseries 
form 

F(z) = a 0 f ( °\z) + a ie - tl2 / (1) (,z) + a 2 z (z) + ..., (2.7) 

where the exponentially suppressed terms 2 correspond to non-perturbative corrections 
with (possibly complex) “instanton actions” ti, while the asymptotic power series f^ l \z) ~ 
£ n>0 /n^ 2 _n_1 correspond to the perturbative expansions on top of the non-perturbative 
corrections 5 . The complex numbers <7* are called transseries parameters, they are constants 
only in wedges of the complex plane, 9\ < arg(z) < 9 2 , but they can jump precisely when 
arg(z) crosses a Stokes line. 

As we will see in a concrete example in Section 4, from the sum over Hankel’s con¬ 
tours as in Figure 1, we can compute the discontinuity of the analytic continuation of the 
perturbative expansion /(°) ( z ) and this is given by 

(<v - Sg-)f(°\z) = A 1 e~ t ' z f( 1 \z) + A 2 e- t2Z fC 2 \z) + ... (2.8) 

where the only ti appearing are the one corresponding to all the singular points in the 
direction 9, i.e. arg ti = 9 , the constants Ai are complex numbers usually called Stokes 
constants (or analytic invariants) while the f^\z) are precisely the asymptotic formal 
power series in 1/z associated to the non-perturbative corrections e~ tiZ . 

In a similar way, the asymptotic series f^‘\z) will have singular directions with non¬ 
trivial Stokes automorphisms that will involve the other expansions z). In this way the 
ambiguities in resumming separately one by one each asymptotic series f^(z) combines 

5 Despite all of the above series being only asymptotic and the transseries expansion being a formal 
coding of the exact function, there is a precise, algorithmic way to compute numerically all the coefficients 
fn once we are given the exact function F(g) with arbitrary numerical precision. We thank Slava Rychkov 
for discussions on this problem. 
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together with the jump in the transseries parameters, and the transseries expansion F(z) 
is an exact representation of the exact answer via Borel-Ecalle resummation 


SoF(z ) = ^ +) S e+ fW(z) + 4 +) e _tl z S e+ f {1 \z) + a { + ] e~ t2Z S e+ f^(z) + ■■■ 

= o-o _) s e- / (0) 0) + C7 1 _) e_<1 z<s e- / (1) (z) + 4 -) e ~ t2 z Se- / (2) («) + ■■■» 


(2.9) 


where the complex numbers are the transseries parameters in the respective wedges 
arg (z) > 8 + e, arg(z) < 0 — e. Even if we do not have at our disposal the full transseries, we 
can still obtain a great deal of information simply by studying the discontinuity properties, 
across a Stokes line, of the resummation of /(°)( z ), or similarly the large orders behavior 
of the perturbative coefficients . 

Let’s suppose that the only Stokes line for f^°\z) is the positive real axis, 0 = 0, then 
by Cauchy theorem 6 [28] we know that 


f W (z) 


1 

27ri 


da; 


/ (0) M 


OJ — z 


_L r d ^. Disco / (0) M 

2vri J 0 uj — z 


( 2 . 10 ) 


where we dropped the contribution coming from infinity, and by expanding for z —> oo 


U3 — Z 


E 

71=0 


Of 1 Z^- 1 


we get 

/n° } ~ — 27ri J 0 nDisc o/M> ( 2 - n ) 

where we schematically wrote f^°\z) ~ 

We can now use (2.8) to relate the large order coefficients /„ to the lower order 
coefficients of the asymptotic expansion f^\z) ~ EE>o fn'z~ n ~ l (we refer to the thorough 
work of Aniceto, Schiappa and Vonk [41] for a more general discussion), then 


f(°) ~ — 
Jn 2vri 


dw u/‘ 


Al e~ tlUJ I H-I + A 2 e- t2 “ I dL_ H-I + • • 


f(i) jy 
Jo i J1 


(i) 


—1 2 w 


/( 


( 2 ) 






a; 


( 2 . 12 ) 


which leads to 

A\ n! 


r(o) 

«/ n 


2vri (ti) n+1 
A? n! 


f (i) , , f d) 

Jo Ji n + J 2 


t 2 


+ /; 


(i) 


l 3 

r i 


+ 


2vri (t 2 ) n+1 


#(2) , XI) *2 , X2) 

JO ' j\ ' j 2 / 1\ 

n n[n — 1) 


n(n — 1) ' J3 n(n — l)(n — 2) 


+ 


t 2 

r i 


(2.13) 


This beautiful result tells us that the large order coefficients of the perturbative expan¬ 
sion do contain explicitly the lower order coefficients of the non-perturbative contributions 
expansion (and vice versa)! Furthermore we want to stress that the large order behavior 

6 Note that generically we would get contributions coming from all the discontinuities across every singular 
directions in the Borel plane. 
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of the perturbative expansion can also be use to determine whether the problem at hand 
is resurgent or not [42]. Had we forgotten to include some non-perturbative contribution 
[43-45], or some transseries parameter (see Section 5), we would have found that the large 
orders of the perturbative expansion do not satisfy (2.13). We can “experimentally” check 
that, to obtain the correct asymptotic behavior for the perturbative coefficients f n , new 
additional transseries parameters and/or new instanton types have to be added to our 
original transseries ansatz! 


3 Solving the BES equation 


The starting point of our analysis is the BES equation [11]. In this section, we review 
how to solve this equation at strong coupling. In particular, we would like to compute the 
strong coupling expansion as high as possible. For this purpose, we follow the approach in 
[22, 35] (see also [46, 47]). 

The BES equation is an integral equation for a rapidity density of a twist-two operator 
with large spin. It is derived from the asymptotic Bethe ansatz equations. The explicit 
form of the BES equation is given by 


a(t) 


t 

e* — 1 


K(2gt,0)-4g 2 



d t! K(2gt, 2 gt')a(t') 


(3.1) 


where the kernel K(t, t') consists of two parts: the “main scattering” part and the “dressing 
phase” part: 

t') = K 0 {t, t') + Ki(t, t'), ^ 

roc y.// 

K d (t, 0 = 8 g 2 / d t" Ki(t, 2gt") K^lgt", t'), 

Jo e - 1 

with 


K (, - t'Jo(t)Jl(t') jv t'Ji(t)J 0 (t') - 

A 0 (M) =--, A!(M) =-- 


(3.3) 


The unknown function a(t) is the density to be solved. Once this function is known, the 
cusp anomalous dimension is given by 


r c us P (5) = 8 c/ 2 cj(0). 


(3.4) 


Therefore our task is to solve the BES equation (3.1). We stress that the BES equation 
is believed to be valid at any coupling. It smoothly interpolates between the weak cou¬ 
pling regime and the strong coupling regime. At weak coupling, it is not hard to solve 
(3.1) directly order by order in g 2 [11],' while to find the strong coupling solution is not 
straightforward. A direct approach to solve the BES equation at strong coupling (and at 
finite coupling numerically) is found in [18, 19]. Below, we review the method developed 
in [22, 35]. This approach is indeed very powerful to find the strong coupling solution 
systematically. 

'However, it is more convenient to use a discrete matrix formulation in subsection 3.2 [18] in the sys¬ 
tematic higher order computation. 
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3.1 Strong coupling solution 

To find the strong coupling solution, we first divide a(t) into even/odd parity parts: 


e*-l ^ 7+(2^) + 7-(2^) 

-a{t) = — 


t ' ' 2gt 

Then the BES equation leads to the following equations: 


7±(-f) = ±7±(*)- 


di [ 7-(2gt) 7+(2gt) 

t 1 — e _t e* — 1 
d t r 7+(2g*) _ 7- (2gt) 
1 — e _t e* — 1 

Next we introduce two functions T±(t) by 


«^2n— l(2gt) — r,^nl; 

'hn{2gt) = o. 


r(t) = r + (t) + ir_(t) = + icoth (7+W + W-W)- 

As was shown in [35], the equations (3.6) is rewritten as the form 

poo 

/ d t [e iuf r_(t) - e~ iut T + {t)} = 2, (|u| < 1). 


(3.5) 


(3.6) 


(3.7) 


(3.8) 


The important point is that the solution T(i) must have an infinite number of zeros at 
t = 47rig(m — 1/4) (m E Z) and poles at t = Airigm' (m! E Z\{0}). This is a reflection 
of the factor 1 + icothin (3.7) and the entireness of <j(t). This analyticity condition 
uniquely fix the solution of the equation (3.8). In general, the equation (3.8) alone does 
not fix its solution uniquely. According to [35], a solution of (3.8) without an analyticity 
condition is generically given by 


P (it) = fo{t)V Q {t) + 


(3.9) 


where Vo(t) and V\ (t) are known functions, whose explicit forms are presented in ap¬ 
pendix A (see (A. 6)). The unknown functions fo(t) and f\{t) are fixed after requiring the 
analyticity condition. The cusp anomalous dimension is then given by 

r c usp(g) = —2gf (0) = 2 5 (1 - 2/!(0)), (3.10) 

where we have used Vo(0) = 1, Vi(0) = 2 and the non-trivial equation (see [35]) 

/o(0) = -1. (3.11) 

The condition that T(t) has the zeros at t = 4irig(m — 1/4) (m E Z) leads to a quantization 
condition 

fo{tm)Vo(t m ) + fi(t m )Vi(t m ) = 0, t m = 4ng fm - ^ . (3.12) 

For later convenience, we rewrite it as 

fo{4TTgx m ) + r(x m )fi(4:TTgx m ) = 0, x m = m - (3.13) 
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Table 1 . First few non-perturbative contributions present in r(x m ) for different values of m £ Z. 


-2 

-1 

m = 0 

+1 

+2 

A 0 

A° 

A° 

A° 

A° 

A 18 

A 10 

A 2 

A 6 

A 14 

A 36 

A 20 

A 4 

A 12 

A 28 


where 


r(x) = 


V\ (Airgx) 


(3.14) 


F 0 (4 ngx)' 

Note that this quantization condition is valid at arbitrary coupling. In [35], the strong 
coupling solution of the quantization condition (3.13) was constructed. One important 
remark is that the asymptotic behavior of r(x) at strong coupling is different for x > 0 and 
x < 0 (and also for Reg > 0 and Re g < 0). This is clearly understood by looking at the 
relation (A.9). It is also important to note that r(x m ) has a non-perturbative correction 
of order 


r np {x m ) ~ 0(Al 8m - 2 '), me Z. 


(3.15) 


See (A.13) and (A.14) for full detail. Therefore, up to 0{ A 4 ), only r np (xo) contributes to 
the quantization condition, while at 0(A 6 ), r np (.xi) also appears. In table 1, we show for 
small values of m, the order in the non-perturbative scale A, for which r(x m ) is going to 
contribute. 

We want to stress in here that the transseries expansion we construct will only be 
valid for Re 5 > 0, simply because the transseries expansion for r(x) has been obtained in 
this half plane, see equations (A.9). Even if the Borel-Ecalle resummation of the complete 
transseries could give us an analytic continuation of the cusp anomalous dimension for 
Re g < 0, this will not be the physical cusp anomaly for Re g < 0. While we know 
from the weak coupling convergent expansion that T cusp (c/) = r cusp (— g), the analytic 
continuation of the transseries would break this symmetry! The correct way to proceed is 
to obtain the transseries expansion for r(x) for Re g < 0 and construct from here the new 
transseries in the left half plane. When we glue together the two Borel-Ecalle resummations 
of the two transseries for Re g ^ 0, we will obtain the physical cusp anomaly satisfying 
Tcusp (. 9) = T C usp(— g), and clearly at the gluing line, the imaginary direction, we will have 
branch cuts, as expected from weak coupling results . 8 

Since r(x m ) receives exponentially suppressed corrections in g —y 00 , the functions 
fo(t) and fi(t) also receive the non-perturbative corrections of the form e~ 2n9 at strong 

8 It is already mentioned in [11, 18] that the r cU sp(<7) has branch cuts along the imaginary axis starting 
at g = ±i/4. 
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coupling. They are thus given by transseries expansions 


m = /o 0) (o + + (£) Vm + °< a6 )> 

2 (3.16) 

AW = /iT) + ^/l (1) (t) + (^) /f’W + 0(A 6 ), 

where the non-perturbative scale is given by (1.6). We want to compute the 1 /g corrections 
in ) (j = 0,1). As was computed up to 0(l/g 2 ) in [35], the perturbative part takes 
the following beautiful form 


where 


/o 0) (4vr gt) = 


e=o 


(4 irgY 




oo 


E 


1 

(4 TrgY 


T><f) p oJ ()) 

70(*)d? ()) 


7o (t) 


r d)r(i -t) 

r(f -t) 


7i(*) = 


+ 7i wqS (j) 

+ 7iWQi°i (j) 

r(j)r(i + t) 
r(| + i) 


(3.17) 


(3.18) 


A remarkable observation 9 is that (1 ft) and Q^J(l/t) (j = 0,1) are polynomials ofl/t 

with degree l. One can check this observation by computing /^(t) and order by 

order, following [35]. Here we assume this observation in all orders in the 1 jg expansion. 
At the leading order: t = 0, we have 



(3.19) 


The polynomials P^(l/t) and (j = 0,1) are uniquely fixed by the following four 

conditions 


• The quantization condition (3.13). 

• The analyticity condition (3.11). 

• The function f\ (47r gt) must be finite at t = 0 because the cusp anomalous dimension 
(3.10) is a finite quantity. 

• The polynomials Q^jil/t), P^j (1 /t) and Q^\(l/t) do not have constant terms, i.e., 

Q$$(o) = A ( ?(o) = Q ( i°l(o) = o. 

Here the last condition is based on the observation in the low order computation. We do 
not have a clear reason of this assumption. We indeed computed the polynomials pj^(l/t) 
and numerically order by order up to n = 180, starting with (3.19). Once these 

9 We thank Benjamin Basso for telling us this remarkable structure. Note that a similar property is also 
found in a different approach in [47]. 
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functions are found, one can immediately compute the strong coupling expansion of the 
cusp anomalous dimension by (3.10). We give the explicit results up to 1 /g 10 in appendix B. 
Since the solution here is based on the above assumptions, we need to confirm its validity. 
We have confirmed the strong coupling expansion of r cusp (g) computed in this way is in 
perfect agreement with the one computed in another way in [47] up to 1 /g 42 . 

In the similar manner, one can fix the non-perturbative corrections /(j "' 1 (t) and f[ n \t) 
in principle. We observe that the leading and the next-to-leading corrections take the 
following form: 


OO 

= Y. 7W 


^=0 

oo 


( 47 rgY 


fl^dngt) = £ 


e=o 


7o(*) p d? (j) +7i(«)Qo?< (j) 

7o(*) p i? (j) +7i(j) 


(n = 1,2), (3.20) 


where 


7o (t) = 7o (t), 


7i (t) 


r(|)r(i + f) 
r(f + t) 


7i (*) 

1 + 4+ 


(3.21) 


The functions p\ n )(l/t) and Q[ n )(l/t) (j = 0,1; n = 1, 2 ) are again polynomials in 1 It with 
degree l. For f = 0we have 


p(4) 

^ 0,0 


P, 


( 2 ) 


0,0 





(?) 

o® (?) 





(3.22) 


We indeed computed p\ n ){l/t) and Q[ n )(1/t) (n = 1, 2) up to i = 180 numerically . 10 These 
solutions will be used to analyze the large order behavior of the cusp anomalous dimension 
in section 5. 

The correction at O (A 6 ) is much more complicated. As mentioned before, at this order 
the quantization condition (3.13) receives the new contribution from r{x\). This additional 
non-perturbative contribution makes the solution more involved. In appendix C, we work 
out the leading correction at 0( A 6 ) by following the original method in [35]. 


3.2 Numerical evaluation at finite coupling 

Here we review how to evaluate the cusp anomalous dimension numerically at finite coupling 
[18]. The basic idea is to reduce the problem from the integral equation to an infinite 
dimensional linear system. The two functions 7 ±(t) in (3.5) admit the following Neumann 
series: 

OO OO 

7-(*) = 2 Y, (2m l)72m —1 +2)71 —1 (+i 7+(+ — 2 ^ ] (2m)72 m ./ 2 m(+1 (3.23) 

m= 1 m= 1 

10 In the non-perturbative sectors, we need to relax the fourth condition above. We require only P( T ° (0) = 
0 for n = 1,2. We confirmed that this condition and the first three conditions above uniquely fix the solution. 
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Table 2. The numerical values of the cusp anomalous dimension. 


9 

r C us P (9)/( 2 <?) 

9 

r c usp(5)/(2g) 

0.1 

0.19385434324817514169 

0.2 

0.35843733036950590918 

0.3 

0.48610560959171757729 

0.4 

0.58210341218853482707 

0.5 

0.65393902847754809138 

0.6 

0.70804410773964657108 

0.7 

0.74930396158229752404 

0.8 

0.78126610325935190716 

0.9 

0.80645604710807850846 

1.0 

0.82665913262694175149 

1.1 

0.84313878800192894948 

1.2 

0.85679498059481188901 

1.3 

0.86827541056713620211 

1.4 

0.87805219572881873579 

1.5 

0.88647437750685446821 

1.6 

0.89380385613927175199 

1.7 

0.90024008980819058115 

1.8 

0.90593719988552639946 

1.9 

0.91101593465281368193 

2.0 

0.91557213314609428817 


Plugging these expansions into the integral equations (3.6), one obtains an infinite number 
of linear equations 

OO 1 

'In T ^ ( fc"nm r )m — ^ 1. (3.24) 

m= 1 


where 




Jo 

Therefore the coefficients 7 n is formally given by 


d t J n (2gt)J m (2gt) 
t e* — 1 


(3.25) 


_ 1 ( 1 
ln 2V1 + /C 


nl 


The cusp anomalous dimension is finally given by 


(3.26) 


r cU sp(sO — &9 Ti 


V 


1 + /C 


11 


(3.27) 


Of course, it is hopeless to compute the inverse matrix (1 + /C ) _1 analytically at arbitrary 
coupling, but in the practical computation it is sufficient to truncate K to a finite dimen¬ 
sional matrix. The convergence is quite rapid as a cut-off of the matrix size grows. For 
example, to get the value at g = 1 with 20-digit precision, it is sufficient to take the 30 x 30 
truncated matrix. In table 2, we show the several numerical values of r cusp (<5r)/(2g) for 
1/10 <g< 2. 


3.3 The mass gap in the 0(6) sigma model 

It is known that the cusp anomalous dimension at strong coupling is closely related to the 
mass gap in the 0(6) sigma model. This was first discussed by Alday and Maldacena in 
the dual string consideration [36] . Then in [33-35] , the relation was embedded into J\T = 4 
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SYM. In particular, in 
BES equation by 


[35], the mass gap mo( 6 ) is exactly related to the solution to the 


m O(6) 


16\/2 /l(-TTff) 

vr 2 Vb(-7rg)’ 


(3.28) 


where /i(t) is the same function appearing in the BES solution. This mass gap scales 
non perturbatively as e -7rs at strong coupling (see (A. 12)). Remarkably, the leading non- 
perturbative correction to the cusp anomalous dimension is completely captured by this 
mass gap [35]: 


r CU s P (sO 

2g 


oo p(0) 


(3.29) 


The mass gap itself also has the following transseries expansion 


m O(6) 


m 


(n) 

0 ( 6 ) 


vi 

r (!) 

OO 

E 


(2vr^) 1/4 ( 


- 1*9 


m 


(o) 

0 ( 6 ) 


(n) 


i=0 


(27 rg) e 


A 2 (i) 

Bv Tg m °(V 


+ 



m 


( 2 ) 

0 ( 6 ) 


+ 0 { A 6 ) , 


(3.30) 


where only the very first few exact values of m!f^ and rVp are found in [35]. One important 
consequence of (3.29) is the following. The “perturbative” coefficients rn f- are computed by 
the perturbative BES solution f[°\t) at t = —ng. Using the relation (3.29), the coefficients 
r{ 1] are fixed by m'p. More explicitly, we have 


°° p(l) 

1 1 

h ( 2 *s)‘ 


E 


L £=0 


m 


(0) 


(27 Tgf 


(3.31) 


We conclude that all the coefficients rjj 1 ' 1 can be predicted only from the perturbative 
solution to the BES equation. As we will see in section 5, we find that a similar relation 
also holds in the next-to-leading non-perturbative sector (see (5.34)). 


4 A toy model 

In [35], Basso and Korchemsky constructed an interesting toy model solution of the integral 
equation (3.8). This toy model solution remarkably captures many features that the actual 
cusp anomalous dimension has. Thus it is useful to understand T cusp (< 7 ) at strong coupling 
from this toy model, as we will see here in detail . 11 

4.1 Transseries solution 

The solution is obtained by relaxing the analyticity conditions (3.12). The toy model 
solution does not have any poles, but has only the zero at t = — 7 ri g. This analyticity 

11 However, one should keep in mind that there are some significant differences between them. We will 
comment on them in the next section. 
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condition can be easily imposed directly in equation (3.9) and one finds 

/<«’”(() =-i, 

/! tor) (t) = -c <,0,l (9), cf'oy'ns) = (4,1) 

Surprisingly, both fQ° y \t ) and /-[ toy ^(f) do not depend on t. The exact solution to the 
BES equation thus takes the form 

r (toy )(ii) = -Vo(t) - c^ toy ) (g)V\(t), (4.2) 


where Vo(t) and V±(t) are the same functions as before, and can be expressed in terms of 
Whittaker functions (see (A. 7)) . 

From the above expression we can obtain the toy model cusp anomalous dimension 



-2 5 f (toy) (0) = 2 g 


1 - 2 


Vo(-Trg) 


2 g 


1 - 


Vi(-n g) 

M 1/4,1/2(27TS0 


V^g m -i/4,o( 27T 9) - ’ 


(4.3) 


where we used the relation (A.7) to the Whittaker function of the first kind. Similarly to 
the true cusp anomalous dimension, the small coupling expansion in the toy model also 
has a finite radius of convergence |<?| < 0.7966..., while the strong coupling expansion is 
only asymptotic as we will shortly prove. 

Using the result (A. 12) derived in the appendix we can rewrite the toy model cusp 
anomaly 

r5(s) 5 + F(i,f|a)+A 2 + 5 0 F(-4,||-a) 1 

2 9 a S + F(l,l\a) + lA^aSoF(ll\-ay ° 2vr ff ’ 

where F(a,b\z) is defined by the asymptotic series (A.l), and S± = 5 0 ± are the lateral 
Borel resummations along the 0 = 0 direction. The Borel resummation of F(a,b\z) is 
defined by (A.2). Similarly, it is also possible to write it as 


r { S(g) = S-F(ll\a) + AlS 0 F(-ll\-a) 

2 5 a 5_F(i,i|a) + tA 2 _a5 0 F(l,H-a)’ 1 j 

where A± are the non-perturbative scales 

A| = TiM'/V 2 '*, (4.6) 

V4/ 


Note that the expressions (4.4) and (4.5) are equivalent, and both Borel resummation 
procedures S± leads to the same result unambiguously thanks to the non-perturbative 
parts A^_. 

We can now expand (4.4) or (4.5) at strong coupling (i.e. for a small with Rea > 0) 
and obtain the transseries expansion for the toy model cusp anomaly 


ituTp ) (ff) 


2 g 


C 0 (a ) - aA 2 C 2 {a) + |a 2 AjU 4 (a) + 0{A \), 
Co (a) - aA 2 C 2 (a) + \a 2 A A _ C A {a) + 0( A 6 _), 


0 <arga<f, 
— \ < arg a < 0 , 


- 15 - 

















where the perturbative expansion Co (a) is 


Cq (a) = 1 — a 




while each function C 2n (a ) is given by 


C 2 (a) 



C 4 (a) 


Fill l«) 3 ’ 


C 2n (a) 


(nf.fl-q ))"" 1 

(F(bl\*)T +1 


(4.8) 


(4.9) 


To obtain these results, we used the Wronskian like relation (A. 5). Given the series ex¬ 
pansion (A.l), it is clear that each C 2n (a), once expanded for small a, is given by an 
asymptotic power series, non-Borel summable for a £ M + . Note however that F (|, ||—a) 
is instead Borel summable for real and positive a while it becomes non-Borel summable 12 
for a £ . 

Following the results derived in Section 2, the Borel-Ecalle resummation of the transseries 
(4.7) is given by 


rS } (g) = f cS+C 0 (a) - aA 2 S + C 2 {a) + \a 2 A\ S+C 4 (a) + • • • , 0 < arga < § , 

2 g 1 S-Co(a) — aA 2 _S-C 2 (a ) + \a 2 A^_S-C^{a) + • • • , — | < arga < 0 . 

(4.10) 

Despite the direction arga = 0 being a Stokes line for both Co and all the C 2n , the 
resummation of the transseries expansion (4.10) is perfectly good and analytic in the entire 
right half-plane — 7 r /2 < arga < 7t/2. As we will show, the ambiguity in the resummation 
prescription (5+ — S-)Cq ^ 0, and similarly (5+ — S-)C 2n ^ 0, is exactly balanced by the 
jump in the transseries parameter cr+ ~ e 3m / 4 —y (J_ e + 37 Tl/ 4 _ 

Despite the fact that the Borel-Ecalle resummation of the transseries (4.10) defines 
an analytic continuation for Fcu°)p (g) in the entire complex plane, minus the negative real 
axis, nonetheless this is not the correct analytic continuation for Re g < 0. Similarly to 
our previous discussion in Section 3, this transseries ansatz is only a solution to the toy 
model quantization conditions in the right half plane Re g > 0, since we used the transseries 
form (A. 12) for Vo an d V\ only valid in this half plane. We should first obtain the correct 
transseries in the left half plane Reg < 0, and then perform its Borel-Ecalle resummation 
to obtain the toy model cusp anomaly. In particular the correct analytic continuation for 
Tcu°p)(g) will not have a branch cut on the negative real axis similar to the true r cusp , as 
one can easily see from the weak coupling analysis. 

Naively, to compute the discontinuity (5+—5-)Co, one expands (4.8) for small a, com¬ 
putes its Borel transform, and then studies the analyticity properties of B[Cq\. However, 
it is much better to use the fact that S± are good resummation prescriptions, i.e., 


5+C 0 (a) 


5+ 


1 — a 


F(bl\a) 
F (l ||a) 


= 1 — a 


S + F( 

S + F( 


i 5 
4’ 4 

1 1 
4’ 4 


a) 

a) 


(4.11) 


12 This means that, while for Co and C 2 the only Stokes line is arga = 0, for higher C 2 n also arga = 7r 
will be a Stokes line. 
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One can indeed check this equality by performing the lateral Borel(-Pade) resummation 
on both hand sides. 

At this point we can make use of the trivial identity 5 + = + (5+ — <S_), together 

with the known discontinuity (A.4) for F(a,b\a), to obtain 

S+C 0 (a) = S_C 0 (a) + 4E (-f^f 

n =1 ' v “ 


S-C 2n (a ), (4.12) 


with Ag = e ±37n / 4 Aj_. The discontinuity of the resummation of the perturbative expansion 
is precisely of the form (2.8), discussed in Section 2. Similar results can be derived for all 
higher perturbative corrections to the non-perturbative terms A^_ n , all schematically of the 
form 

* A 2 

S+C 2n (a) = S-C 2n (a) — (fl + 1)-^= S-C 2 ( n +l)(oi) + ^(Aq) (4-13) 

We can now understand why the Borel-Ecalle resummation (4.10), despite looking discon¬ 
tinuous for arg a = 0, defines indeed an analytic function in the right half-plane: the change 
in the transseries parameter in front of the C 2 term is taking care of the non-perturbative 
contribution of order Ag in the discontinuity (4.12) of Cq. The jump of the C 4 coefficients 
balances the sum of the non-perturbative Ag term in the discontinuity of Co together with 
the non-perturbative contribution of order Ag coming from the discontinuity of C 2 (4.13), 
and so on. 

Furthermore, from the discontinuity equation (4.12) for Co (a), we can use the result 
(2.13) for the large order coefficients 13 of the perturbative expansion to derive 


C , 


2 A 


°’ (n) ^vri(l)™- 1 / 2 
2 A 2 


F n ~ 2 


1 + 


a 


2 ,( 1 ) 


X 1 


n~ | 


+ ^2,(2) x l 2 , 0 (n - 3 } j + 

(n-§)(n- |) I 


7ri (2) ? 


—j-T (n — 1) I 1 + 


C 4 ci'i x 2 C 4 ( 2 ) x 2 2 , , \ 

] -xr + 0(n -3 ) + (4.14) 


n — 2 


{n — 2){n — 3) 


+ 


2 A 3 


7 ri (3) n 3 / 2 


r In -- 1 + 


c, 


6 ,( 1 ) 


x 3 


n~ I 


+ 0(n 


— 2 \ 


+ • • • 5 


where C 0 ( n ), C 2j t n ), ■■■ are the n-th order perturbative coefficients of the small a expansion 
of Cq, C 2 ,..., while the Stokes constant A is simply given by 


A = 


—ir(|) ilmcr + — iImcT_ 


2 V / 2T(|) 


(4.15) 


as one could have also deduce from the jump in the transseries parameter (see [39]). 

As mentioned in Section 2, the large order behaviour (4.14) of the perturbative coeffi¬ 
cients C 0i ( n ) can be used to test the validity of our transseries expansion. From equation 
(4.8), we can easily generate an arbitrary number of perturbative coefficients C 0j ( n ), and 

13 Unlike the case in the previous section, there is a fractional exponent factor a in A§. This factor 
shifts the arguments of the Gamma function, but the essential argument does not change at all. 
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check the leading “one instanton” 14 behaviour C 0 ^ ~ Ar(n — 1/2). We can proceed 
by considering the first 1/n correction to the leading T(n — 1/2) term, this will give us 
the first perturbative coefficient, C- 2 . m, of the one instanton sector, that can be checked 
against the analytic value obtained from the expansion of (4.9). Higher and higher terms 
C' 2j ( 2 ), C 2j ( 3 ),... can be obtained by simply looking at higher n~ 2 ,n~ 3 ,... corrections. 

We already know that the C 2 expansion is also asymptotic, so the coefficients C 2 ( n ) 
will diverge as well. For this reason we can perform a Borel-Pade resummation of the 
leading one instanton sector, llnst ~ AT{n — 1/2)(1 + C 2 i m/(n — 1/2) + ...), and subtract 
it from the perturbative coefficient C 0 r n y In this way we can isolate the subleading large 
order behavior (4.14) which is C(w re ) — llnst ~ A 2 Y{n — l)/2 n . From this new set of 
coefficients we can read the perturbative coefficients C 4i of the two-instantons sector, 
repeat the Borel-Pade resummation and isolate the three-instantons sector and so on. 15 

We checked numerically, up to the four-instantons sector, that the large orders of 
the perturbative series expansion do indeed contain all the non-perturbative information, 
consistently with our one parameter transseries expansion (4.7). 


4.2 The non-perturbative mass scale 


In Section 3, we briefly discussed the relation between the non-perturbative scale of the 
cusp anomalous dimension and the mass gap of 0(6) sigma model [36]. A similar analysis 
[35] can be carried out for the toy model as well. In the toy model the non-perturbative 
mass is exactly given by 

16\/2 1 . . 

mtoy ~ 7 r 2 Vi (- 7 rg)' (416) 

Using the results derived in the Appendix (A. 12), we can rewrite the mass gap as 

mtoy= 7rT(|) a 7 6 2 “5 ± T(|,i|a) + |aA^ 0 F(|,||-a)' (4 ‘ 17) 

Note that in the toy model the mass gap takes a particularly simple transseries form, 

we just need to expand equation (4.17) around A± = 0, and we obtain a one parameter 

transseries. The only reason why we get infinitely many non-perturbative corrections in 
rntoy is that we expand 1/Ui(— ng), had we considered l/mt oy ~ Vi(— ng) we would have 
obtained a very simple two-ternrs transseries (see [39]). 

As in (4.7), we expand the mass gap as 


^toy — 




-a 




a A 2 a 2 At 

m t 0 y,o{a )- oy ,i(a) H-^^rn toy , 2 (a 0 + ©(A^) 


(4.18) 


for 0 < arga; < 
given by 


For 0 < arga < |, one replaces A + —> A_. 


^■to y,n (®) 


[F(ll\-a)] n 


(n> 0). 


The coefficient mt oy ,n is 


(4.19) 


14 In this paper, we sometimes use the term “instanton” to count the non-perturbative order A 2 for 
notational simplicity. Keep in mind that this does not mean that such corrections are caused by instantons 
in the 0(6) sigma model. There is no instanton configuration in this model, and the Borel singularities 
correspond to IR renormalons. 

15 We thank Marcel Vonk for useful discussions on this problem. 
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It is easy to see that there are non-trivial relations between C 2 n {oi) and mt 0 y,n, 

^ltoy,n(oO ^toy,l(oO -^( 4 , 4 ! 

m toy,n—l{ a ) m toyfi( a ) (4.20) 

C 2 n(cx') — ?7ltoy,0 (oO^ttoy.n—1 (^0 • 

The latter relation with n = 1 is the same as the one in (3.31). In the toy model, the 
similar relations to (3.31) hold for all the non-perturbative sectors. In the next section, we 
also confirm that the similar relation hold for n = 2 in the true cusp anomalous dimension. 

5 The cusp anomalous dimension at strong coupling 

In this section, we study the strong coupling expansion of the physical cusp anomalous 
dimension. We first see singularities of the Borel transforms. Then we analyze the large 
order behavior of the perturbative expansion in detail. The result shows that the infor¬ 
mation on the non-perturbative sectors are encoded in the perturbative sector. We next 
perform the lateral Borel resummation, and compare it with numerical values of T cusp (g) 
computed directly from the BES equation. 

5.1 Borel singularities 

Let us first see the singularity structure of the Borel transform. The cusp anomalous 
dimension has the transseries expansion (1.5). Using the method reviewed in section 3, 
we have computed the numerical values of rj ”' 1 (n = 0 , 1 , 2 ) up to £ = 180 with 200 -digit 
precision. A very first few values are in perfect agreement with the exact ones in [22, 35]. 
Some higher order corrections are presented in appendix B. 

Since T^ grows as T(^ — 1 / 2 ) in l —> oo, it is natural to consider the following a bit 
modified Borel transform 

OO -p(ra) 

BFSpKC) := E fvr-r^' P- 1 ) 

Since we have only the finite number of coefficients T'/’’ up to l = 180, we need an ap¬ 
proximation of S[Tcusp](C)- A natural way is to replace it by its Pade approximant. In the 
following, we use the diagonal Pade approximant of £>[r£usp](C) with order 90. 

In figure 2, we show the Borel singularities of the Pade approximant for £>[riui p ](£) 
(n = 0,1,2). The figure clearly shows that the Borel transforms I?[Pcusp](C) f° r n = 0,1 
have the singularities at ( = 1, —4, while £>[Pcusp](C) has the singularities at £ = ±1. Note 
that in a Pade approximant, a condensation of poles indicates a branch cut of the original 
function. Therefore it is very likely that £>[Pcusp](C) for ra = 0,1 has two branch cuts 
(—oo,—4) and (l,oo) while B[r/ u ; sp ](C) has branch cuts (—oo, — 1) and (l,oo). 

These singularity structures are very important to understand the form of the full 
transseries, as emphasized in section 2. The large order behavior of T), ' also heavily 
depends on the Borel singularities. 

Note that from figure 2 it is manifest that also the direction 9 = n in the complex 
Borel plane is a Stokes direction. Following the discussion of Section 2 it is clear that 
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Figure 2. We show the singularities of the Pade approximant for the Borel transform f?[r£usp](C) 
(n = 0,1,2) in the complex Borel plane. It is obvious to see that the Borel transforms for n = 0,1 
have the singularities at £ = 1, —4, while the Borel transform for n = 2 has the singularities at 

C = ±1. 


from the large order behavior of the perturbative coefficients we will see the presence of 
these discontinuities. However by the point of view of the cusp anomalous dimension, this 
additional Stokes line is not relevant. In fact as we have already reviewed in Section 3, the 
transseries that we are studying, is only valid in the right half plane. It is not necessary to 
add a transseries parameter to take into account the new singularity at £ = —4, and the 
Stokes line 6 = n because it is outside of the regime in which the Borel-Ecalle resunmration 
of our transseries reproduces the physical cusp anomaly. Of course the large order behaviour 
of the perturbative coefficients will know nonetheless of this additional singularity in the 
negative axis of the Borel plane. 

5.2 Large order behavior 

We will study now the large order behavior of the perturbative coefficients. First of all, 
as was observed in [22], the perturbative coefficients rj 0 '* grow factorially as T(l — 1/2) in 
the large i limit. As in section 2, the large order behavior of rj, 0 '* has more information. It 
contains the non-perturbative information. In fact, following the argument in section 2, it 
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is expected that T 
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behaves in the large l limit as 16 
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-2 (£-2)(£-3) 
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+ ■■■ 


oo. 


(5.2) 


where the first two coefficients Ai and A 2 are exactly the same as the toy model ones in 
(4.14)-(4.15): 


2 A 

Ai = — = -0.3042971194498708318670259057666.. 

7T1 

2 A 2 

A 2 = -= 0.14545061420433549651130121616 ... i. 

7ri 


(5.3) 


As we will see later, the constant A 3 is more subtle. Below we want to confirm the relation 
(5.2) by using the numerical values of T^°\ T^ and Pp computed from the BES solution. 
We first estimate A\ . Let us consider a sequence 


,(!) _ 


40) 




(5.4) 


This sequence must converge to A\ in i 
shows 

a (1) r 
u 180 


00 . The naive evaluation of aP at £ = 180 


-0.304049. 


(5.5) 


This is indeed close to the exact value in (5.3), but the agreement is not so good because of 
the sub-leading contribution. In fact, the sub-leading contribution is of order l/£ ~ 5x 10 -3 . 
To improve the precision, we need to remove such corrections. This can be done by using 
the Richardson extrapolation. Let us consider the ra-th Richardson transform [48] of a 
sequence fg, 

(- l) k+n (£ + k) r 


Knlfl] ■= E 


fc =0 


k\{n-k)\ 


-ft+k- 


(5.6) 


Assuming (5.2), the sequence a~P behaves as 


af ] = Ai [1 + 0{r 1 )} , l ->■ 00 . (5.7) 

The Richardson transform of aP then behaves as 

n n [a { p] = A 1 [l + 0(r n - 1 )}, 00 . ( 5 . 8 ) 

16 As shown in figure 2, the Borel transform f3[riui p ](^) has the singularity at £ = —4. This singularity 
implies a contribution of order (— 4)~ e to the large order behavior. Therefore this contribution is important 
to see r) 4) from the large order behavior. 


- 21 


















Difference 



Figure 3. (Left) We show the convergence behavior of the original sequence af \ its first and 
fifth Richardson transforms by the blue, purple and red curves, respectively. The black dashed line 
is the expected convergence value A\. The Richardson extrapolation accelerates the convergence. 
(Right) The difference 11 — 1Z n [a^go_ -rJ/^il * s s h° wn as a function of n. In this case, the Richardson 
extrapolation is especially good for the range 20 < n < 60. 


Therefore the convergence speed is improved. Note that to compute the n-th Richardson 
transform of a ( 1> we need the higher coefficients ,..., o!fl n - Thus if we have up 
to certain hnite l = £ max , we can perform the n-th Richardson transform only up to 
£' = l ma.Y — n. It is a trade-off how we should choose n and £'. 

In the left of figure 3, we sketch the convergence behavior of the Richardson extrapola¬ 
tion for the sequence a*p ■ We plot the original a!f \ its first and fifth Richardson transforms 
up to £ = 20 by the blue, purple and red solid curves, respectively. It is clear to see that 
the Richardson transform indeed accelerates the convergence. We also show, in the right 
figure, how the Richardson transform lZ n with fixed £ max = 180 works as n grows. In the 
current case, one should choose n in the range 20 < n < 60. In the following analysis, we 
use the 30th Richardson transform. Applying the Richardson transform to the sequence 
a) , we get 

^ 30 [ag 0 ] « -0.304297119449870831867025905573, (5.9) 


which is in agreement with (5.3) with 27-digit accuracy! 

Once A\ is fixed, we can proceed to the next-to-leading coefficients T^. We consider 
a sequence 




do) 


AT(^-l) 


- 1 


(5.10) 


This should converge to T) 1 in £ —> oo. After the Richardson transform again, we find 


7^30[&i5o] ~ -0.144860385419958982062924086. 


The exact value of T 


(i) 

l 


is found in [35]: 


(5.11) 



3(1 - 2 log 2) 


-0.144860385419958982062924091 


8 


(5.12) 













Table 3. The numerical estimation of from the large order behavior (LOB) in the perturbative 
expansion. We use up to i = 180 and the 30th Richardson extrapolation. We also show the 
deviation from the true r^ 1 -*. 


~i Estimation from LOB |T^-7^3o[—j/T^ 


1 

-0.1448603854199589820629240863 

3.3 x 10 -26 

2 

-0.0938694159959333439731081 

5.5 x 10 -23 

3 

-0.14973248501597359266824 

5.3 x 10 -21 

4 

-0.402728449526514898734 

2.2 x 10 -19 

5 

-1.4747555976090392157 

5.9 x 10 -18 

6 

-6.86483498151787623 

1.0 x 10" 16 

7 

-38.674657843660018 

1.3 X 10" 15 

8 

-255.9115123079947 

9.8 X 10" 15 

9 

-1945.413350672867 

1.7 x 10" 14 

10 

-16710.570082722 

9.3 x 10" 13 


Both precisely coincide as expected. Pushing the same computation, we can systematically 
confirm the relation (5.2) order by order. In table 3, we summarize the results estimated by 
the Richardson extrapolation. The large order behavior indeed captures the information 
on the leading non-perturbative sector. 

Next, we want to test the second line in (5.2). To do so, we have to subtract the 
contribution in the first line. This can be done by using the Borel resununation. We first 
rewrite the first line in (5.2) as 


p(i) 

' !<I T> : = 1 + 

Z 2 






i + E 

m= 1 



(5.13) 


The important point is that the coefficient 7^ grows factorially, and thus h^(£) is a 
formal divergent series. Furthermore, we observe that the Borel transform t3[h^](C) has 
a singularity at £ = log 2. Thus is non-Borel sunnnable, and one has to use the 

lateral Borel resummation. The ambiguity of the two lateral Borel resununations is roughly 
estimated as 

(5+ - 5_)/r (1) (£) ~ i x 0(2~ e ), (5.14) 

where S± = 5 0 ±. other words, the imaginary parts of S±h^\£) are of order 0(2~ e ). 
Now we consider a quantity 

drf := rj 0) - Air (e - 0 s+h^(£). (5.15) 

As mentioned above, since the imaginary parts of S±h^(£) behave as G(2~*), dT / (l ' i behaves 
as i x 0(2~ : ) in £ —> oo. This just corresponds to the contribution in the second line of 
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Table 4. The numerical estimation of 1’^ from the large order behavior of r^°\ We use I^ 0 ' 1 up 
to t = 180 and the 20th Richardson extrapolation. 


£ 

Estimation from LOB 

i -77.20 [■■•]/rf ) | 

1 

-0.7499999999999987 

1.7 x 10" 15 

2 

0.845472324053663 

7.0 x 10" 14 

3 

-2.4877562953846 

8.8 x 10" 13 

4 

7.504881881284 

7.2 x 10" 12 

5 

-37.76786450152 

1.1 x 10" 11 

6 

186.98364073 

8.6 x 10" 10 


(5. 2). 17 Let us consider a sequence 



r(*-i)' 


(5.16) 


This should converge to A As before, using the Richardson extrapolation, we find 


7*20[agi] ~ 0.1454506142043354998i + 3.0 x 10“ 2 °, (5.17) 


whose imaginary part is in agreement with the exact A 2 up to 17-digit accuracy. As in 

( 2 ) 

the leading non-perturbative sector, one can estimate the coefficients T) ' in (5.2). The 
obtained values up to £ = 6 are shown in table 4. 

Finally, let us proceed to the third line in (5.2). At this level, it is not easy to extract 
the information with high precision by numerics. We again consider quantities 


2 T 


( 2 ) 


* (2) W :=l + ^ + 


2 2 r 


( 2 ) 

2 


-2 (£ — 2)(£ — 3) 


+ 


= 1 + E 


( 2 ) 
7m 


m =1 


(5.18) 


and 

d 2 r{ 0) := rf - Air (e- 0 s+hP-\e) - A 2 r ^~ 1 1) s + h^(£). 

Defining two sequences 18 


we find 


(3),odd _ 




T(2£- 1- f 


_a 2 t (0) 

" 1 2£—1 ’ 


(3),even _ 


5 2 i-i 


T(2£- f) 


(5 2 r (0) 

0 1 2 1 i 


77 8 [a£ ),odd ] w 0.101563618725, 
77 8 [4o ),even ] w 0.101563618713, 


(5.19) 


(5.20) 


(5.21) 


17 The relation (5.2) can be regarded as a “transseries” in the large l expansion. The second line is 
the leading “non-perturbative” correction in 1 /£, whose “instanton action” is log 2. The ambiguity of the 
lateral Borel resummations of the perturbative part (the first line) is almost canceled by the ambiguity in 
the second line. In other words, we can guess the contribution in the second line from that in the first line. 

18 The reason of the separation of the sequence is purely technical. We observed that the Richardson 
extrapolation works better for this rather than for the naive sequence. 
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Unlike A\ and A 2 , these do not agree with the Stokes coefficient in the toy model: 

2A 3 

— = 0.0695237641771486153 .... (5.22) 

7T1 

This discrepancy is explained as follows. As noted in section 3, in the quantization condition 
(3.13) up to order A 4 , only the non-perturbative correction of r(x m ) with m = 0 gives a 
non-vanishing contribution. However, at the order A 6 , r np (x m ) with m = 1 also contributes 
to the quantization condition. As we will see in appendix C, if taking into account this 
contribution, we find the exact value of * 4.3 as 

„ 2A 3 ( 8T 4 (|) \ 

„4 3 =- 1 + = % = 0.101563618709385381..., 5.23 

vri V 3V3T 4 (f )J 

which is indeed very close to the numerical estimations (5.21) from the large order behavior! 
This is precisely the expected value of the Stokes constant to reproduce the jump of the 
two transseries parameter responsible for the A 6 corrections, see the Appendix C, equation 
(C.12). This result implies that at the order A 10 , the Stokes coefficient also receives a new 
contribution from r np (x m ) with m = —1. 

The reason for this different behavior that begins only at the three instantons level 
0(A 6 ) can be understood from the quantization condition. The quantization condition 
(3.13) is a set of infinitely many linear equations in the undetermined functions fo(t) and 
fi(t). The coefficients appearing in these equations are the known functions r(x m ) and at 
strong coupling they assume the transseries forms (A.15)-(A.16). If we were to turn off all 
the non-perturbative corrections in all the r(x m ), then our strong coupling perturbative 
expansion f^°\t), in (3.16), would solve this infinite set of equations. 

We can proceed to turn on the first non-perturbative correction in the quantization 
condition (3.13). As shown in table 1, the first non-perturbative correction appears at 
order 0( A 2 ) in r(x o). This means that a non-perturbative correction of the same order 
has to be added to and f^\t) and this correction has to be constructed entirely 

from the particular expansion (A.15)-(A.16) combined with the perturbative expansions 
as in [35]. The next non-perturbative correction is of order 0( A 4 ) and it is 
once again coming from the transseries expansion of r(x’o). A correction of the same order 
will now be needed in and /W(i) and this correction has to be constructed from 

the order 0( A 4 ) of r(x 0 ) combined with our perturbative expansions and 

from the order 0( A 2 ) of r(x 0 ) combined with the 0( A 2 ) just obtained for 

From table 1, we see that at order 0( A 6 ) we will receive two contributions, see (A.15)- 
(A.16), one from r(x 0 ) and one from r(x 1 ). The correction of this order to and 

f^\t) can still be obtained from the lower orders, but it stops being a one parameter 
transseries and becomes a two parameter transseries, see Appendix C. Clearly, as we see 
from table 1 , at 0( A 10 ) the non-perturbative corrections to r(x 1 ) will come into play and 
we will need yet another transseries parameter and so on so forth. 

The above large order behavior, as well as our analysis of the quantization condi¬ 
tion equations, strongly indicates that the strong coupling expansion of the physical cusp 
anomalous dimension is a multi-parameter transseries, while the expansion in the toy model 
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is a one-parameter transseries. This is a big difference between the toy model and the true 
cusp anomalous dimension. 


Large order behavior of the mass gap. Let us turn to the large order behavior of 
the mass gap coefficients m\ . Using the results derived in the Appendix C, we computed 
the large order behaviors of the perturbative expansion in (3.30), as well as for the 

coefficients of the first two instantons corrections and . 

0 ( 6 ) 0 ( 6 ) 

As expected from the discussion in Section 2, the perturbative coefficients m^ behave 
as 




+ 


b r(<-)) 

2jri 3'- 1 / 2 


(i + o(r 1 )) + ..., 


a 2 r(£-i) 
+ 2 £ ~ 1 



2m 


( 2 ) 


1-2 


+ 0{£ 


- 2 \ 


(5.24) 


where the Stokes constant A is precisely the one obtained in the toy model (4.15) or in the 
cusp anomalous dimension. 

The Stokes constant for the three instantons sector B is once again not consistent with 
a one parameter transseries ansatz, i.e. it is not simply B ^ A 3 , see also the discussion at 
the end of Section 5.2. At the three instantons level a new transseries parameter arise and, 
from the correct multi parameter transseries form derived in the Appendix (C.13), we can 
deduce 

d = 17r 

8\/3r(f) 2 ' 

For the large order coefficients in the first non-perturbative correction in (3.30), 

we can easily recognize the presence of the perturbative coefficients as well as the higher 
instantons one: 


m 


(i) 


c r(e + l) / , —m 


(o) 


27ri 


1 + 


l - - 
c 2 


+ 0 {£ 





£-3/2 


+ 0 {£ 


—2\ 


(5.26) 


where the Stokes constant C is associated with the presence of a Stokes line in the 6 = n 
direction, with a characteristic alternating factor 1 /(—2 tt) £ in front of it, and it is given by 


C = 



(5.27) 


The presence of this Stokes line is obvious since the one instanton sector will “see” the 
perturbative sector with a relative action 27t x (—1), hence the alternating nature of the 
coefficients. Note that this is precisely the same Stokes constant in the 9 = n direction 
that one can obtain from the toy model! 
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5.3 Lateral Borel resummations and ambiguity cancellation 

In this subsection, we perform the Borel resummation of the strong coupling expansion. As 
emphasized in [22], the perturbative strong coupling expansion of r cusp (g) is a non-Borel 
summable asymptotic expansion. This is because that the Borel transform of Tcusp (. 9 ) bas 
the singularities on the positive real axis in the Borel plane, as shown in figure 2. One 
has to avoid these singularities when performing the inverse Borel transform (or Borel 
resummation). There is an ambiguity to choose such integration contours, and the lateral 
Borel resummations are not real-valued any more. The imaginary part must be canceled by 
the other contributions, i.e., the non-perturbative sectors, in the full transseries expansion 
(1.5) because the cusp anomalous dimension is a real-valued quantity for g £ M. Here we 
test this cancellation by explicitly performing the lateral Borel resummations of Tc”sp(g) 
for n = 0,1, 2. 

Let us consider the inverse Borel transform of (5.1). By using the well-known integral 
representation of the gamma function: 

roo 

T(z)= dCe" c C 2 "\ Rez > 0, (5.28) 

Jo 

it is easy to see that the inverse Borel transform is given by 

s 0 rW p ( 5 ) = 1 + J™ d(e-<c 3 / 2 B[r£> p } (J-'j 

roc ~ 

= 1 + (2vn?)- 1 / 2 / dCe- 2 ^C- 3/2 S[r^ p ](C) 

Jo 

where we have separated the l = 0 term from the others to use the integral representation 
(5.28). This resummation is, however, not well-defined because H[Fcusp](C) has singularities 
on the positive real axis. Therefore we need to modify it as the lateral Borel resummations: 


(5.29) 


S±r <&(</) = 1 + ( 2 t rg )- 1 / 2 


dCe -2^C c -3/2£ [r (n) ](C) 


(5.30) 


These lateral Borel resummations are complex-valued quantities, and cause the disconti¬ 
nuity (5+ - 5_)riui p (5f). 

Now we perform the lateral Borel-Pade resummations at finite coupling. Let us see in 
detail for g = 1/5 as an example. Using the formula (5.30), we find the following numerical 
values 

5 ± r iu Sp (l/5) ~ 0.1252604 0.2097335i, 

5 ± rW p (l/5) « 0.8678150 =f 0.08111881, (5.31) 

5 ± r iu Sp (l/5) ~ 0.6406029 =F 0.0846409i. 

These are complex numbers, and obviously the imaginary part in the perturbative part is 
not negligible. More importantly, these resummations take the different values depending 
on the integration contours. Therefore there is an ambiguity of the resununations. Plugging 
these values into the transseries (1.5) and setting the parameter as a = e^ 371 " 1 / 4 [35] for 
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Table 5. The ambiguity cancellation in the lateral Borel resummations at g = 1/2. The total 
imaginary part is indeed very close to zero due to the cancellation. The real part should be 
compared with the directly evaluated values by the method in section 3. 


Contributions Values 

5±r^ p 0.63182936 =f 0.02188690i 

- ft 5 ± r ^sp 0.02210593 ± 0.02166913i 

j^pS ± T ( H p 0.00000326 ± 0.00021822i 

Sum 0.65393855 ± 0.00000044i 

Direct Evaluation 0.65393903 


Table 6. The ambiguity cancellation at g = 1. 


Contributions Values 

S ± riu } sp 0.825965638690 =f 0.000693266968i 

- H 5 ± r cus P 0.000693493818 ± 0.000693040133i 

I g^S±r£ ) sp 0.000000000111 ± 0.000000226842i 

Sum 0.826659132619 ± 0.000000000007i 

Direct Evaluation 0.826659132627 


respectively, we finally get 

•Sitfip - ||s ± rWp + 

« (0.1252604 T 0.2097335i) + (0.2303182 ± 0.1909411i) + (0.00249307 ± 0.01886876i) 

» 0.3580716 ± 0.0000763i, (5.32) 

where A^_ are defined by (4.6). The imaginary part becomes very small after taking into 
account the non-perturbative corrections! Furthermore, the real part is close to the true 
value directly evaluated from the BES equation in table 2. We also show the ambiguity can¬ 
cellation for g = 1/2,1, 2 in tables 5, 6 and 7. The real parts are in remarkable agreement 
with the values in table 2. We conclude that the ambiguity of the lateral Borel resumma¬ 
tions are precisely canceled by the perturbative and non-perturbative contributions and 
that the finial results are unambiguous real values. 

5.4 Novel relations 

In this subsection, we observe quite novel relations between the cusp anomalous dimension 
and the mass gap in the 0(6) sigma model. So far, we do not have a derivation of these 
relations from the BES equation. It would be important to prove them along the line [35]. 
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Table 7. The ambiguity cancellation at g = 2. 


Contributions Values 

«S±r^ p 0.915571204186344585360 =p 0.000000928959331963843i 

- ^5 ± r^ S p 0.000000928959749702527 ± 0.000000928958914225168i 

rcusp 0.000000000000000000282 ± 0.000000000000417738679i 

Sum 0.915572133146094288169 ± 0.000000000000000000005i 

Direct Evaluation 0.915572133146094288174 


We first find that the strong coupling expansion of the mass gap (3.30) has the following 
beautiful relation 


m O(6) 

(i) 

m O(6) 


m, 


(i) 

0 ( 6 ) 


m. 


(o) 

0 ( 6 ) 


0.6051396146 0.8516809822 

1---+ 


2ng 


(2 ngf 


2.271453077 
(2vr g) 3 + 


(5.33) 


We checked this equality up to 1/g 180 with sufficient numerical precision. Also, we find the 
similar relation to (3.31): 


oo 


E 



oo (0) 

\ - rn\ 

' oo (1) 

\ - 

' (2ngY 
u =o v J 

2^ (2TrgY 
u=o v J 


(5.34) 


These relations suggest us the following guess: 


m, 


(n) 

0 ( 6 ) 


(n-1) 


m 


0 ( 6 ) 


(1) 

? m O(6) 

OO 

V 

p( n ) ? 

oo (0) 

oo (n— 1)" 

(0) ’ 
m O(6) 

e=o 

( 2-KgY 

2^ l2ngY 
li=Q y J 

2^ (2 ngY 
Lf=o v J 


(5.35) 


In the toy model in the previous section, we have already seen that the similar relations 
hold for all n. However, in the physical case, this is probably not the case for n = 3. 
A reason of this is very likely due to the contribution from r np (xi). 19 This contribution 
should change the structure of the coefficients T) , but the right hand side in (5.35) there is 
no source of such a change. It would be interesting to check it in more detail by computing 
the higher non-perturbative corrections from the BES equation. 


6 Conclusions and future directions 

The cusp anomalous dimension has a rich non-perturbative structure at strong coupling. 
We can put it on the resurgent analysis, recently developed in many contexts. Solving 
the BES equation at strong coupling systematically, we computed the perturbative and 
non-perturbative corrections up to very high orders in the 1/g expansions. Using these 

19 In other words, if we turn off this correction by hand: r np (a:i) = 0, then the relations (5.35) seem to 
work for n = 3. Of course, this does not give the true cusp anomalous dimension. 
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data, we showed that the large order behavior in the perturbative expansion has the non- 
perturbative information, as expected. Moreover, the ambiguity of the Borel resumma¬ 
tions in the perturbative sector is precisely canceled by the contributions from the non- 
perturbative sectors. The final result is real-valued and in remarkable agreement with the 
direct evaluation at finite (not so strong) coupling. We also find the unexpected relations 
between the cusp anomalous dimension and the mass gap of the 0(6) sigma model. 

In this paper, we focused on the cusp anomalous dimension. Clearly, it would be 
interesting to explore the strong coupling expansions for other examples. Closely related 
examples are the generalized scaling function, proposed in [49], and the generalized cusp 
anomalous dimension (or equivalently the quark-antiquark potential), studied in [50-54], 
As studied in [33, 55-57], the strong coupling analysis of the generalized scaling function 
is almost in parallel with the cusp anomalous dimension, and thus it is a good exercise to 
see its resurgent aspect along the line in this paper. The analysis of the generalized cusp 
anomalous dimension will be much more involved. 

Another example is the eigenvalue of the adjoint BFKL operator. This eigenvalue plays 
a very important role in scattering amplitudes in the so-called multi-Regge limit. Recently, 
Basso, Caron-Huot and Sever proposed integral equations that compute the adjoint BFKL 
eigenvalue at any coupling [58] (see also [59]). Since these equations are quite similar to 
the BES equation, it seems to be possible to use the technique in section 3. It would be 
interesting to ask whether the BFKL eigenvalue receives non-perturbative corrections at 
strong coupling or not. 

The strong coupling limit of short operators is also an important problem. In partic¬ 
ular, the Konishi operator is a significant example. In [60, 61], it was observed from the 
TBA analysis [62-65] that the conformal dimension of the Konishi operator has the strong 
coupling expansion 


A k - 2 = 


( 


1 H— 1= + • • • 
\/A 


oo. 


( 6 . 1 ) 


The same result was also reproduced from the various worldsheet computations [66-69]. 
It would be important to clarify the analytic structure of the Konishi operator at strong 
coupling. The Quantum Spectral Curve formulation [5, 6] will be helpful for this purpose. 
See [70] for an interesting approach to the strong coupling Konishi dimension, based on 

[71]- 

The generalization to the ABJM theory [13] is also interesting. The all-loop Bethe 
ansatz equation in planar ABJM theory was conjectured by Gromov and Vieira in [72] . One 
important consequence is that in the integrability approach, the’t Hooft coupling always 
appears through a non-trivial function h{ A), called the interpolating function. Then, it 
turned out in [72] that the cusp anomalous dimension in ABJM theory is related to that 
in AT = 4 SYM by 

r^ M (A) = ^= 4 (<? = K A)). (6.2) 

It was a long-standing open problem to determine h{ A). Surprisingly, Gromov and Sizov 
recently proposed an exact formula of h( A) [15]. According to their result, the interpolating 
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function at strong coupling behaves as 


MA) = ~ ~ir + °( e ' 27Tv ^)’ ( 6 - 3 ) 

Now it is clear to see that T^?p M (A) receives two kinds of non-perturbative corrections at 
strong coupling. One is just obtained by replacing g in r^= 4 (g) by h{ A): 

0 ( e - 2 nh W) ~ 0{e~ nV ^). (6.4) 

The other is the correction 20 that h{ A) itself receives, as in (6.3), 

0(e~ 2nV ^). (6.5) 

It would be interesting to discuss physical origins of these two corrections in the effective 
worldsheet theory on AdS 4 x CP 3 in the high spin limit. As discussed in [14], the latter 
type correction is probably related to stringy worldsheet instantons, constructed in [75]. 
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A Special functions 

Let us introduce some functions used in the main text. See appendix D in [35] in more 
detail. We first consider the following series: 


n«.&l*):=EwW^‘- (A-l) 

k =0 ^ ' k 

where we used the Pochhammer symbol (x) n = T(x + n)/T(x). This series is asymptotic. 
Let us consider the following “Borel transform”: 

^(o = E7if|[ = ( 1 -cr 6 - (a.2) 

k =0 ^ ,k 

This Borel transform has a branch cut along (l,oo) if b is a non-integer. Therefore for 
z < 0, F(a,b\z) is Borel summable. Since the series is non-Borel summable for z > 0, one 

2(, This non-perturbative correction is quite different from the first type correction (6.4) even though 
the exponential factors are almost same. The strong coupling perturbative expansion of h(X) in (6.3) is 
a convergent series. Thus one does not need the Borel resummation. Nevertheless, h( A) receives non- 
perturbative corrections! This non-perturbative structure is essentially the same as that in the 1/2 BPS 
Wilson loop [73], analyzed in [14, 74]. In fact, the idea in [15] is to relate the interpolating function to the 
results in [14, 74]. 
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needs the lateral Borel resummations. One easily finds that the inverse Borel transform of 
(A.2) is given by 


I i-ooe ^ 10 

S 0 ±F(a,b\z) = —— / dCe _c C a_1 (l - z()~ b . 

r (a) Jo 


(A.3) 


A simple computation shows that the discontinuity of these lateral Borel resummations is 
given by 


V 1 —a—b 


(S 0 + - S 0 -)F(a, b\z) = 2vri- 


_1 / 2 5oF(l — a, 1 — b\—z), z > 0. 


mm 

Remarkably, the Borel resumnrations satisfy the following Wronskian-like relation [35] 


(A.4) 


S «* F l i j 




_z I “ I i' i 




~Z = 1, 


(A.5) 


where we assumed z > 0. 

Next, let us define the functions V n (x ) and U^(x) (n = 0,1) by 


V n (x) 

Ut{x) 


— ^ du(l + u) 1 / 4 -"(l -um 4 e ux , 

k J -1 


1 


d u (u ± l)~ 1/4 (u l) 1 / 4-n e _ ( u_ (Re x > 0). 


(A.6) 


Note that the integral in V n (x) is well-defined for all complex values of x, but the integrals 
in U±(x) are convergent only for Rex > 0. For Rex < 0, we need analytic continuations 
of Un(x). As mentioned in [35], the function V n (x) is written in terms of the Whittaker 
function of the first kind 

UM = ^ r(l)r(2- B ) ( ~ 2:c) " /2 " lMl/4 -" /2 ' 1/2 -^ 2( ~ 2;,:) ' (A ' 7) 

Also, the functions U±(x) are related to the Whittaker function of the second kind: 


u m = (|) ^ 1 e"VF_ 1/ 4 , 1/2 ( 2 x), U+(x) = Q) ( 2 x)- 1 / 2 e a; W 1/ 4 ,o( 2 x), 

M*) = (i) x- 1 e"W 1/4 , 1/2 (2x), ur(x) = Q) (2x)- 1 / 2 e"W_ 1/4 ,o(2x). 

(A.8) 

These relations naturally give analytic continuations for Re x < 0. Note that these functions 
have branch cuts on the negative axis. As shown in [35], the functions V n (x) and U^{x) 
are related by 

V 0 {x) = ^e^mUoi-x) +e~ x U+(x)}, 

2% (A ‘ 9) 
V\(x) = —e T37ri/4 [e a: C/ 1 -(-x) - e~ x U+(x)\, 

7r 


where the upper and lower signs correspond to Imx > 0 and Imx < 0, respectively. For 
real x, one needs the ie-prescription on the right hand sides due to the branch cut of U^{x). 
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Now we rewrite U^(x) in terms of the Borel resummation of F(a,b\z). As in [35], for 
x > 0 we have 


K(x) 

Uo(x) 

U+(x) 

ur(x) 


= (2x)- i l i T (!) S 0 F (j, ||- 
= (2x)^‘r(?)*F(-i,» 

= i(2 I )- 1/ Jr (i) 50 F(i,i 



(x > 0), 


(A.10) 


As mentioned before, F(a,b\z) is Borel summable for z < 0. If the argument of U±(x) is 
negative, we need the ie-presciption. It just corresponds to the lateral Borel resummations 
of F(a, b\z) for z > 0: 


K (-* ± ie ) 

Uq (-X ± ie) 
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= ^e =F3 ^ i/4 (2x)- 3 / 4 r ( 


1 

2x 

1 
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V 4 


(x > 0), (A.11) 


Substituting these expressions into (A.9), one can express V n (x) in terms of F(a,b\z). In 
particular, one hnds 
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(A.12) 


where the non-perturbative scales A± are defined by (1.6) with a = a± = e : F 37rl / 4 ) respec¬ 
tively. Note that these expressions for upper and lower signs are equivalent because of the 
discontinuity (A.4) of the lateral Borel resummations. In other words, the ambiguity of the 
lateral Borel resummations is precisely canceled by the non-perturbative correction (the 
second term in (A.12)). 

It is useful to write down the expression (3.14) with x = x m = m— \ at strong coupling. 
As noted before, the strong coupling behavior is different for m> 1 and m < 0. For m > 1, 
we obtain 

1 F(l ili^r) - (-l) m 4(4m - 1)5 (g||) 4m a 2™-l A 8m-2 F (l ; 1]-^) 

r(Xm) = 2- ! /r^n 4m -’ 

H-h + (-l) m (4 m - l)-3 (^j a 2m A 8m-2 F (l, ||__^_) 

(A.13) 
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where a = 7 ^. Similarly, we find the expression for m < 0, 


r(^) = 


2(1 - 4m) 
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^(l, !It 4^) + - 4m)-§ (£[fj) ^ a _ 2 m+ 1 A 2 _ 8 m F(|, fh^) 


^(l, lira) + - 4m) i (§|) ^ a- 2 ™A 2 - 8 -F(-i, f h^) 

(A.14) 

It is obvious to see that the non-perturbative correction to r{x m ) starts from 0(Al 8 m_ 2 l). 

With the above expression (A.14) for r(x m ) we can easily obtain its transseries expan¬ 
sion. Let us assume that Rea > 0, then for m > 0 we can write 
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(A.15) 


where the choice on transseries parameter a± depends on Im a > 0 or Im a < 0 (note that 
for a = l/(27rg) this corresponds to Im g < 0, or Img > 0 respectively). Similarly for 
—m < 0, assuming always Rea > 0, we have 


r(x- m ) = 


2(4 m + 


b \^4 ’ 4 1 (4m+iy J n=l V V ' >J 
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(A.i 6 ) 


where once again the choice of transseries parameter a± inside A± is correlated with Im a > 
0 or Ima < 0 . 


B Strong coupling expansions 

Here we explicitly give the strong coupling expansion up to the next-to-leading non- 
perturbative sector . 21 The perturbative strong coupling expansion is 
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1 More precise and higher order coefficients are available on request to the authors. 
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These are, of course, in agreement with the result in [22]. The expansion around the leading 
non-perturbative sector is also given by 
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Similarly, the “perturbative” part of the mass gap in (3.30) is 
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The expansions in the non-perturbative sectors are given by 
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(B.2) 


(B.3) 


0.07243019271 0.04955777441 0.07845572166 0.2082747743 


(B.4) 


(B.5) 


(B.6) 


C Three instantons corrections 


As we mentioned in the main text, the solution to the quantization condition (3.13), and 
hence the cusp anomalous dimension, stop being a one parameter transseries at the order 
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0(A 6 ). To show this we analyzed the solution to (3.13) using the methods of [35]. As 
shown in [35], we can rewrite (3.13) as 
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with meZ, where the undetermined functions fo(t) and fi(t) are written in terms of the 
coefficients c±(n,g) 
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These coefficients can be found by using the transseries ansatz 
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and solving order by order in 1 /g and order by order in A the infinite set of equations 
(C.l). 

Then, proceeding in the same way as [35], we find the first few coefficients. In the 
perturbative sector, the solutions are 
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In the leading non-perturbative sector, we find 
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At the next-to-leading order, we get 
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Up to this order all of the above coefficients are perfectly consistent with a one pa¬ 
rameter transseries expansion for the cusp anomalous dimension, essentially because up to 
this order we could have truncated to the perturbative level all but the m = 0 equations 
in (C.l). At the three instanton level though, the transseries expansion (A.15) for r(x i) 
contains a A 6 term. 
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Using all of these results, we get the following strong coupling expansion of r cusp (g): 
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These values are all consistent with the results presented in Appendix B. 

/Q\ 

Note that the leading order of r'c U ; sp (g) is different from one, this means that the order 
0 ( A 6 ) should actually be rewritten as 
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with the usual transseries parameter a± = e =F37rl/,4 r(|)/r(|) and the new second transseries 
parameter 


a ± = - e ^ 3 ^ 1 / 4 
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3V3T(f) 

From the correct transseries expansion for the coefficients c±(n, g), we can easily obtain 
the mass gap transseries expansion from equation (3.28). As it turns out also the mass gap 
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is described by an transseries with infinitely many transseries parameters and at order A 6 
is given by 
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Note precisely how the A 6 leading order is not the one that could be guessed from a one 
parameter transseries ansatz! 

Obviously also these values are all consistent with the results presented in Appendix 
B. 
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